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Chaos, randomization, and turbulence in particle-laden flows 


A. Bershadskii 
ICAR, P.O. Box 31155, Jerusalem 91000, Israel 


The randomization effect of the two-way (particle-flow) interaction has been studied and quanti- 
fied using the notion of distributed chaos and the results of numerical simulations and laboratory 
measurements. It is shown, in particular, that an increase of such parameters as the particle volume 
fraction, particle mass loading, and Stokes number results generally in stronger randomization of the 
particle-laden flows. An important role of spontaneous breaking of the local reflectional symmetry 
in the randomization of the particle-laden flows has been also analyzed using relevant dynamical 


invariants. 


I. INTRODUCTION 


A general classification of the non-laminar regimes in 
fluid dynamics can be based on the concept of smooth- 
ness. Smooth dynamics is characterized by the stretched 
exponential spectra 


E(k) x exp —(k/kg)?*. (1) 


where 1 > 6 > 0 and & is wavenumer. The particular 
value 8 = 1 - the exponential spectrum: 


E(k) « exp(—k/k.), (2) 


is usually associated with deterministic chaos (see, for 
instance, Refs. [1]]-[5}). 

For 1 > £ the chaotic-like dynamics is still smooth 
but not deterministic (the distributed chaos). The 
non-smooth dynamics is characterized by the power-law 
(scaling) spectra. 


In the distributed chaos the value of the parameter 
6 can be used as a measure of randomization: a lesser 
value of @ (i.e. it is further from the deterministic value 
8 = 1) corresponds to stronger randomization of the flow. 


In the particle-laden flows the particles can be an im- 
portant additional source of the fluid dynamics random- 
ization (whereas the increase in the Reynolds number is 
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FIG. 1: Kinetic energy spectra computed in the isotropic de- 
cay at t= 5. The spectra are vertically shifted for clarity. 
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FIG. 2: Kinetic energy spectra computed in the isotropic de- 
cay at t = 4 for different values of the particle volume fraction 
gy . The spectra are vertically shifted for clarity. 


still the main source of randomization). Figure 1, for in- 
stance, shows the influence of the different types of parti- 
cles on the fluid’s kinetic energy spectrum. The spectral 
data were taken from Fig. 3 of the Ref. (6). In this pa- 
per results of direct numerical simulation (DNS) of the 
decay of incompressible isotropic flow (with and without 
inertial particles) were reported (we will return to the 
description of this DNS in more detail below). Now one 
can see that for the unladen (particle-free) case the ki- 
netic energy spectrum is exponential (it is indicated by 
the dashed curve), i.e. corresponds to the deterministic 
chaos (the Taylor-scale Reynolds number Re, = 31 for 
the time of decay when this spectrum was computed). 
For the particle-laden cases, the spectra can be well fit- 
ted by the stretched exponentials Eq. (1) (the dashed 
curves), i.e. these cases correspond to the (smooth) dis- 
tributed chaos. The values of the parameter ( are differ- 
ent for the different types of particles and a clear trend 
for randomization can be observed. 

The analogous trend can be observed in Fig. 2 where 
results of an analogous DNS but with different values of 
the particle volume fraction ¢, have been shown (the 
spectral data were taken from Fig. 10 of Ref. [7). 


In the next sections we will relate the values of 3 ob- 
served in the Figs. 1 and 2 with invariants of the fluid 


dynamics and will show an important role of the sponta- 
neous breaking of the local reflectional symmetry in the 
randomization of the particle-laden flows. 


II. THE FLUID DYNAMICS INVARIANTS AND 
DISTRIBUTED CHAOS 


When the characteristic scale k. in the spectrum Eq. 
(2) randomly fluctuates one should use an ensemble av- 
eraging to obtain the average spectrum 


E(k) « 7, P(k.) exp —(k/k.)dke (3) 


If the randomized dynamics is still smooth, then E(k) 
has the stretched exponential form Eq. (1). Compar- 
ing Eq. (1) and Eq. (3) one can find asymptote of the 
probability distribution P(k,) for large k. (1) 


P(ke) 0¢ kz t+8/PO-P) exp(—yhB/C-8)), (4) 
here y is a constant. 


The ideal (nondissipative) fluid dynamics has two fun- 
damental invariants: energy and helicity. The dissipative 
(Navier-Stokes) dynamics also has two fundamental in- 
variants: Birkhoff-Saffman integral {3} {9}, {10} 


Ing = f (u(x,t) wbx +1,1))ar, (5) 
and Loitsyanskii integral Aik | 


Ip= [Plas t)-u(x+r,t))dr (6) 


where < ... > denotes a global (ensemble or spatial) 
average, u(x,t) is the velocity field. 

Due to the Noether’s theorem the conservation of 
the Birkhoff-Saffman and Loitsyanskii integrals by the 
Navier-Stokes dynamics is a consequence of the space 
homogeneity and isotropy respectively. 


Let us denote an invariant as J and relate the charac- 
teristic velocity u. and the characteristic scale k. using 
dimensional considerations 


Ue X 1° Ke (7) 


where 6 and a are some parameters that can be found 
from the dimensional considerations. 


For the normally distributed characteristic velocity ue 
the distribution P(k-) can be readily found from the 
relationship Eq. (7). Comparing this distribution with 
the asymptote Eq. (4) one obtains a relationship between 
the exponents a and 8 


B= (8) 


log, E(k) a.u. 
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FIG. 3: Spectrum for the velocity fluctuations experimentally 
measured at the location behind the grid where the aerosol 
particles, droplets, and the chaotic fluid flow actively inter- 
acted. 


For instance, using the Birkhoff-Saffman integral Eq. 
(5) as invariant J in Eq. (7) we obtain 6 = 1/2 and 
a = 3/2, and consequently from the Eq. (8) 


E(k) x exp —(k/kg)?/4 (9) 


For the flow dominated by the Loitsyanskii integral Eq. 
(6) 6 = 1/2 and a = 5/2, and consequently from the Eq. 
(8) 

E(k) « exp —(k/kg)°/®. (10) 

One can recognize these spectra in the Figs. 1 and 2. 

Let us now discuss the DNS reported in the Refs. fal [j 
in more detail. 

The dispersed solid particles (80 million) were injected 


in the freely decaying isotropic and homogeneous flow 
described by the Navier-Stokes equation 


Ou _ 


1 
a —(u- V)u— Bre Naar Ep (11) 


for an incompressible fluid V-u = 0. The diameter of 
particles was smaller than the Kolmogorov scale. 


The term F, in the Navier-Stokes equation (11) is 


F, = a sh (12) 


here f,, is the drag force acting on particle n and mf is 
the fluid mass within the control volume. 


The particle equation of motion was taken in the form 


dv, u—Vp 


"Ply = Mp Tp + (Mp = my)g (13) 


here v, and u are the instantaneous velocities of the 
particle and of the carrier flow at the particle location, 


d/dt, is the time derivative following the moving parti- 
cle, m, is the particle mass, 7, is the particle response 
time, g is the gravitational acceleration. This equation is 
expected to be a good approximation for heavy particles 


(Pp/p > 1). 


The chaotic motion of the carrier flow was generated 
by initial conditions in the form of a random Gaussian 
noise with an energy spectrum 


Eo(k) x kexp —(k/ko) (14) 


The injection of the particles was done at t = 1 in the 
DNS time terms. 


For the spectral data shown in Fig. 1 the ‘micropar- 
ticles’ correspond to the case of small Stokes number 
Sty = T/T «1 (the 7, is the Kolmogorov time scale), 
the ‘critical’ particles correspond to St, ~ 1, and the 
‘large’ particles correspond to large Stokes number 
St, >> 1 (at fixed particle volume fraction ¢, = 107% 
and particle mass loading ¢,, = 1), the gravity is not 
taken into account. For the spectral data shown in Fig. 
2 the Stokes number was fixed St, = 1 and the particles 
volume fraction , was varying. 


It was observed in the Ref. (6) that the microparticles 
were not centrifugally ejected from the vortical struc- 
tures of the carrier fluid due to their fast response to the 
fluid velocity fluctuations. Therefore, the conditions of 
quasi-isotropy and quasi-homogeneity can be still applied 
to this case. The critical particles were ejected from the 
vorticity cores but were not completely ejected from the 
original vortical structures due to the centrifugal effect 
as the large particles do. The spectral consequences of 
these differences can be seen in Figs. 1 and 2 (cf Eqs. 
(9), (10), and (20)). 


Let us now discuss the results of a recent laboratory 
experiment reported in Ref. (13). In this experiment, 
a moist-air wind tunnel flow behind a passive grid 
with the aerosol particles injected in the measuring 
section was used in order to experimentally simulate 
‘aerosol-cloud-turbulence’ interactions. Since the 
Re, = 30 the flow was not turbulent but chaotic. This 
value of Re, is very close to that corresponding to Fig. 
1 (see above). Therefore, one can compare the results 
of this experiment with those shown in the Fig. 1, 
especially because the flows behind grids are usually 
used to simulate isotropic and homogeneous situations 
(the droplets in the flow were small). 


Figure 3 shows the spectrum for the velocity fluc- 
tuations measured at the location behind the grid 
where the aerosol particles, droplets, and the chaotic 
fluid flow actively interact. The dashed curve indi- 
cates the stretched exponential spectrum Eq. (10) (cf 
Fig. 1) and the dotted arrow indicates the position of kg. 


Ill. SPONTANEOUS BREAKING OF LOCAL 
REFLECTIONAL SYMMETRY 


The previous consideration was concentrated on 
isotropic, homogeneous flows with global (net) reflec- 
tional symmetry. The global reflection symmetry results 
in zero global (mean) helicity. The point-wise helicity, 
however, may not identically be equal to zero in this 
case. Spontaneous breaking of the local reflectional 
symmetry (and related spontaneous helicity fluctua- 
tions) can be considered as an intrinsic property of 
the chaotic/turbulent flows (see, for instance, Refs. 
(14) ,[15),[16]). The emergence of the moving with the 
fluid well-defined vorticity blobs with non-zero_blob’s 
helicity can accompany this phenomenon f17}-[23h. 
Eventually, the vorticity blobs with high relative helicity 
and low viscous dissipation can be a generic property 
of the chaotic/turbulent flows (18). The low viscous 
dissipation in the blobs can result in adiabatic invariance 


of their helicity (17), 20}. 


The well-known phenomenon of inertial particle 
ejection from the areas of high vorticity (strengthened 
by the high relative helicity) can preserve this adia- 
batic invariance for the particle-laden flows. For the 
microparticle-laden flows, in which the particles are 
mainly moving with the fluid (see above), this adiabatic 
invariance should be also naturally preserved. 


Since the global/net helicity should be zero at the 
spontaneous breaking of local reflectional symmetry the 
localized positive and negative blobs’s helicities should 
be canceled at the overall average. 


The helicity in a vorticity blob is 
H; = i h(x, ¢) dr. (15) 
Vj 


where A(r,t) = u-w is helicity distribution, V; is 
spatial volume of the j-blob, u is velocity field, and 
w(r,t) = V x u(r,t) is vorticity field. 


The moments of the helicity distribution h(r, t) = u-w 
can be then defined as ; 


’ 1 rn 
In = Jim 5 ys Hi (16) 
j 
here V is the total volume of the blobs. 


Let H; be the helicity of the blobs with negative he- 


licity and H 7 be the helicity of the blobs with positive 
helicity. Then corresponding sign-defined moments are 


i 
TT 


V—-0o 


* 1 xoyn 
= lim 7 HF (17) 


The summation in Eq. (17) is made over the blobs with 
negative (or positive) Hj only. 


Due to the 
L=tF+Ir = 


global/net — reflectional 
0 for the odd n. 


symmetry 
Then for the 


It should be noted that the blobs with high relative 
helicity provide the main contribution to the moments 
with high n (though, often the value n = 2 can be 
already considered as sufficiently high for strongly inter- 
mittent chaotic/turbulent flows). This gives additional 
support for the use of the helicity moments I, and I* as 


n 
relevant adiabatic invariants for the particle-laden flows. 


IV. DISTRIBUTED CHAOS AND THE 
HELICITY MOMENTS 


Let us begin, for simplicity, with the flows dominated 
by the third moment of the helicity distribution I} as a 
quasi-invariant. 

If one uses the estimate Eq. (7), then from the dimen- 
sional considerations 


te 0 [Eg |S ke!” (18) 


and for the normal distribution of u, one obtains 
P(ke) 


P(ke) oc kz" exp —(ke/ 4h) (19) 
Substitution of the Eq. (19) into Eq. (3) gives 
E(k) ox exp —(k/ke)'/? (20) 


One can compare Eq. (20) with the Figs. 1 and 2. 


In a general case of a flow dominated by the quasi- 
invariant J,, for an even moment or by the quasi-invariant 
I> for an odd moment one can write using the dimen- 
sional considerations 


digo [Tee ee (21) 
for the odd moments and 
Ue  [A/2” Kon, (22) 
for the even moments. Here 
3 
n=l1l-—, 23 
a om (23) 
Then from the Eqs. (8) and (23) one obtains 
2n-—3 
n 24 
ea (24) 


Let us consider two end cases: 
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FIG. 4: Three-dimensional kinetic energy spectra computed 
in the homogeneous free decay at t = 3 for different values of 
the Stokes number St,,. The spectra are vertically shifted for 
clarity. 
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FIG. 5: Kinetic energy spectra computed in a statistically 
stationary isotropic homogeneous fluid motion for Re, = 35.4, 
Stn = 5 (the spectra are vertically shifted for clarity). 


forn > 1 


E(k) x exp —(k/kg)?/3, (25) 


and for n = 2 (i.e. for the Levich-Tsinober invariant 


[19}) 


E(k) x exp—(k/kp)'/* (26) 

Figure 4 shows the three-dimensional kinetic energy 
spectra computed in freely decaying homogeneous fluid 
motion for Re, = 130. The spectral data were taken 
from Fig. 3b of the Ref. (29). In this DNS the gravi- 
tational acceleration g in the Eq. (13) was taken into 
account and the particle-laden flow was anisotropic. The 
Stokes number St, = T/T, is the varying parameter, 
éu =3 x 10-*, and p,/p= 5000. 


The dashed curves indicate the stretched exponentials 
Eq. (10) (for the particle-free case), and Eqs. (25) and 
(20) (for the particle-laden cases). As for the previously 
considered cases of the isotropic free decaying situation 
without gravitation, the situation with gravitation 
exhibits the trend of the randomization enhanced by the 


Ss 

c 

Z| 

WW —$,=0, B=2/3 

> 0) — 9 =0.2, B=1/2 

2 4. — 6,=0.5, B=1/3 
— do =1, B=1/3 


0 1 
10919 k 


FIG. 6: Kinetic energy spectra computed in a statistically 


stationary isotropic homogeneous fluid motion for Re, = 62 
(the spectra are vertically shifted for clarity). 
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FIG. 7: Kinetic energy spectra computed in a statistically 
stationary isotropic homogeneous fluid motion for Re, = 140 
(the spectra are vertically shifted for clarity). 


particles. Unlike the free decaying cases with low Re) 
(shown in the Figs. 1 and 2) in the case shown in the 
Fig. 4 the original particle-free flow was already in the 
state of distributed chaos with the stretched exponential 
spectrum Eq. (10) because the increase_of Re, results 
in enchanted randomization on its own [23}). 


Figure 5 shows the kinetic energy spectra computed 
in statistically stationary isotropic homogeneous fluid 
motion for Re, = 35.4. The spectral data were taken 
from Fig. 9 of the Ref. (24). It should be noted that 
the spectra shown in Figs. 1-4 were computed for the 
free decaying chaotic flows, whereas the statistically 
stationary state was reached in the DNS reported in 
the Ref. using an external random forcing. For 
the particle-laden cases the Stokes number was fixed 
St, = 5, and the particle volume fraction ¢, was varying. 


The dashed curves indicate the stretched exponentials 
Eq. (9) (for the particle-free case), and Eqs. (20) 
and (26) (for the particle-laden cases). As for the free 
decaying situation, the statistically stationary situation 
exhibits the trend of the randomization enhanced by the 
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FIG. 8: Kinetic energy spectra computed in a statistically 
stationary shear homogeneous flow for Re, = 80 (the spectra 
are vertically shifted for clarity). 


particles. Unlike the free decaying cases (shown in the 
Figs. 1 and 2) in the randomly forced case shown in the 
Fig. 5 the original particle-free flow was already in the 
state of distributed chaos with the stretched exponential 
spectrum Eq. (9). Therefore, one can observe in the Fig. 
5 the spectrum Eq. (26) corresponding to the highly 
randomized distributed chaos. 


Figure 6 shows the kinetic energy spectra computed 
in a statistically stationary isotropic homogeneous fluid 
motion for Re, = 62. The spectral data were taken from 
Fig. 5a of the Ref. 25}. The particle mass loadings 
gm is the varying parameter in this case. The dashed 
curves indicate the stretched exponentials Eq. (25) 
(for the particle-free case), and Eqs. (20) and (26) (for 
the particle-laden cases). And again in the randomly 
forced case shown in the Fig. 6 the original particle-free 
flow was already in the state of distributed chaos with 
the stretched exponential spectrum Eq. (25) (there 
was an increase in Re, compared with the previous 
example). Therefore, one can observe in the Fig. 6 
the spectrum Eq. (26) corresponding to the highly 
randomized distributed chaos. 


Figure 7 shows the kinetic energy spectra computed 
in a statistically stationary isotropic homogeneous fluid 
motion for Re, = 140. The spectral data were taken 
from Fig. 2 of the Ref. The Stokes number 
St, is the varying parameter (particle volume fraction 
¢@y = 2 x 10~+, the particle/fluid density ratio is equal 
to 1000). The dashed curves indicate the stretched 
exponentials Eq. (25) for the particle-free case and Eq. 
(20) for the particle-laden cases. 


Finally, let us consider the results of the numerical 
simulations of a particle-laden shear flow reported in 
a recent paper Ref. 27}. A new Exact Regularized 
Point Particle method was applied in the Ref. toa 
homogeneous shear flow with the shear strength S* = 7. 


Figure 8 shows the kinetic energy spectrum computed 


at Ry = 80, St, = 1 and different values of the mass 
loading parameter ¢,, = 0, 0.4, and 0.8. The spectral 
data were taken from Fig. 4 of the Ref. (27. The dashed 
curves indicate the stretched exponentials Eq. (25) for 
the particle-free case, and Eqs. (20) and (26) for the 
particle-laden cases with @m = 0.4 and dm = 0.8 respec- 
tively. 


V. CONCLUSIONS 


As one can see from the above-given examples an 
increase in the particle volume fraction, particle mass 
loading, and Stokes number results generally in stronger 
randomization of the particle-laden flows. The inter- 
play of these parameters provides a rather complex 
picture (value of Re, and presence or absence of the 
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external random force also should be taken into account). 


Relevant dynamical invariants (both dissipative and 
non dissipative) dominate the randomization process and 
the spontaneous breaking of local reflectional symmetry 
plays an important role in this process providing a 
higher level of the randomization. 


The notion of distributed chaos can be instrumental in 
quantifying this phenomenon and the value of the param- 
eter 3 can be used as a measure of the randomization. 
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